Exact results for spin dynamics and fract ionization in the Kitaev Model 
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We present certain exact analytical results for dynamical spin correlation functions in the Kitaev 
Model. It is the first result of its kind in non-trivial quantum spin models. The result is also 
novel: in spite of presence of gapless propagating Majorana fermion excitations, dynamical two 
spin correlation functions are identically zero beyond nearest neighbor separation. This shows 
existence of a gapless but short range spin liquid. An unusual, all energy scale fractionization of a 
spin -flip quanta, into two infinitely massive vr-fluxes and a dynamical Majorana fermion, is shown 
to occur. As the Kitaev Model exemplifies topological quantum computation, our result presents 
new insights into qubit dynamics and generation of topological excitations. 
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In the field of quantum computers and quantum com- 
munications, practical realizations of qubits that are ro- 
bust and escape decoherence is a foremost challenge [ij. 
In this context Kitaev proposed Q certain emergent topo- 
logical excitations in strongly correlated quantum many 
body systems as robust qubits. In a fault tolerant quan- 
tum computation schemejl, 0, 0|i Kitaev constructed 
a non-trivial and exactly solvable 2-dimensional spin 
modelQ] and illustrated basic ideas. In some limit it also 
becomes the celebrated 'toric code' Hamiltonian. The 
Kitaev model has come closer to reality, after recent pro- 
posals for experimental realizations and schemes for 
manipulation and detection^]. In initialisation, error 
correction and read out operations, it is 'spins' rather 
than emergent topological degrees of freedom that are 
directly accessed from outside. Thus an understanding 
of dynamic spin correlations is of paramount importance. 

We present certain exact analytical results for time 
dependent spin correlation functions, in arbitrary eigen- 
states of the Kitaev Model. Our results are non-trivial 
and novel, with possible implications for new quantum 
computational schemes. Further our result is unique in 
the sense that it is the first exact result for equilibrium 
dynamical spin correlation functions in a non trivial 2D 
quantum spin model. 

We show that dynamical two spin correlation functions 
are short ranged and vanish identically beyond nearest 
neighbor sites for all time t, for all values of the coupling 
constants Jx, Jy and Jz, even in the domain of J's where 
the model is gapless. Our result shows rigorously that it 
is a short range quantum spin liquid and long range spin 
order is absent. We obtain a compact form for the time 
dependence, which makes the physics transparent. 

Kitaev Model is known to support dynamical Majo- 
rana fermions and static 7r-flux eigen-excitations. We 
show how fractionization [1, Q of a local spin-flip quanta 
into a bound pair of static tt-Aux excitations and a free 
Majorana fermion occurs. 

In the present paper we have restricted our calcula- 
tion to dynamical correlation functions for time indepen- 
dent Hamiltonians, in arbitrary eigen-states and thermal 



states. In actual quantum computations, key manipula- 
tions such as braiding involve parametric change of the 
Hamiltonian and adiabatic transport of topological de- 
grees of freedomQ. In principle, some of the needed 'non 
equilibrium' dynamical correlation functions may be ob- 
tained by convolution of our results with suitable Berry 
phase factors. 

In our work we follow KitaevQ and use the Majo- 
rana fermion representation of spin-half operators and 
an enlarged Hilbert space. What is remarkable is that, 
because of the presence of certain local conserved quan- 
tities in the Kitaev Model, Hilbert space enlargement 
only produces 'gauge copies', without altering the en- 
ergy spectrum. This luxury is absent for standard 2D 
Heisenberg models when studied using enlarged fermionic 
Hilbert space 

The Kitaev Hamiltonian is 

(*i>x (:ii)y (y% 

where i, j label the sites of a hexagonal lattice, («j)a, a = 
X, y, z denotes the nearest neighbor bonds in the a'th di- 
rection. The model has no continuous global spin symme- 
try. All bond interactions are Ising like, albeit in different 
quantisation directions x, y and z, in three different bond 
types, making the model quantum mechanical. Further, 
it renders a high degree of frustration; that is, even at 
a classical level a given spin can not satisfy conflicting 
demands, from 3 neighbors, of orientations in mutually 
orthogonal directions. The model has a rich local sym- 
metry. A specific product of 6 spin components in every 
elementary hexagon, cr(o-2^^^\'^b'^t (figure 1) commutes 
with the full Hamiltonian. Thus there is one conserved 
Z2 charge ±1, at every dual lattice site of the hexago- 
nal lattice. The model is exactly solvable and becomes 
non interacting Majorana fermions, propagating in the 
background of static Z2 gauge fields. Different possible 
Z2 charges separate the Hilbert space into super selected 
sectors. The ground state corresponds to all Z2 charges — 
+ 1. In this sector, for a range of J's, Majorana fermions 
are gapless, including the special point Jx = Jy ^ Jz- 
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FIG. 1: Elementary hexagon and 'bond fermion' construction. 
A spin is replaced with 4 majorana fermions {c,c^,c^,c^). 
Bond fermion X(23) for the bond joining site 2 and site 3 is 
shown . Spin operators are also defined. 

Following Kitaev, we represent the spins in terms of 
Majorana fermions. At each site, we define 4 Majorana 
fermions, c", a = 0, a;, y, z, 

{c",c^} = 2^„/3 (2) 

Four Majorana (real) fermions make two complex 
fermions, making the Hilbert space 4 dimensional. No- 
tionally, Hilbert space dimension of a Majorana fermion 
is V^) an irrational number, reminding us that Majorana 
fermions have to occur in pairs (leading to a a/2 x \/2 = 2 
dimensional Fock space) in physical problems. 

The dimension of Hilbert space of N spins is 2^. The 
enlarged Hilbert space has a dimension 4^ ~ (\/2 x \/2 x 
y/2 X v^)^. State vectors of the physical Hilbert space 
satisfy the condition, 

A|*)phys = |*)phy. (3) 

^ Q 444 (4) 
The spin operators can then be represented by, 

a°;^icic°;, a = x,y,z (5) 

When projected into the physical Hilbert space, the op- 
erators defined above satisfy the algebra of spin 1/2 op- 
erators, [af,a'j] = iCabcCr'^Sij . The Hamiltonian written 
in terms of the Majorana fermions is, 

a=x,y,z (ij)^ 

with = ic^Cj. Kitaev showed that [iJ, u^y^^] — 

and ui^ij^a become constants of motion with cigcn- values 
~ The variables are identified with static 
(Ising) Z2 gauge fields on the bonds. Kitaev Hamiltonian 
(equation 6) has a local Z2 gauge invariance in the ex- 
tended Hilbert space. For practical purposes, the local 
Z2 gauge transformation amounts to 1*^^)^ TiU(^ij)aTj, 
with Ti ± 1. Equation ([3]) is the Gauss law and the phys- 
ical subspace is the gauge invariant sector. 



In the gauge field sector we have gauge invariant Z2 
vortex charges ±1 (0 and 7r-fluxes), defined as product of 
U(^ij'ja around each elementary hexagonal plaqauette. 

Equation 6, with conserved U{ij)^ is the Hamiltonian 
of free Majorana fermions in the background of frozen 
Z2 vortices or 7r-fluxes. Since Z2 gauge fields have no 
dynamics, all eigenstates can be written as products of a 
state in the 22^ dimensional Fock space of the Ci Majo- 
rana fermions and the (2)'^^ dimensional space of Z2 link 
variables. We will refer to the former as matter sector and 
the latter as gauge field sector. Gauge copies (eigen-states 
with same energy eigen-values) spanning corresponding 
extended Hilbert space are obtained by local gauge trans- 
formations nU(^ij)aTj. 

It turns out that if we attempt to calculate spin-spin 
correlation functions with the use of above free Majorana 
Hamiltonian and the Z2 fields u^^^^'s, it is difficult to 
proceed further. 

It is here we have invented a simple but key trans- 
formation that facilitates exact computation of all spin 
correlation functions. We call this as 'bond fermion' for- 
mation. In the process we also discover a 'quantum frac- 
tionization' phenomenon in the Kitaev Model, that has 
an unusual validity at all energy scales. 

Hereinafter, we follow the convention that i in the bond 
{ij)a, belongs to A and j to B sub-lattice. We define 
complex fermions on each link as, 

xfe% = \ {4 + ^4) (7) 
4). = ^ - '4) (8) 

The link variables are related to the number operator of 
these fermions, = ic°c° = 2x|y)^X(y>„ - 1- All 

eigenstates can therefore be chosen to have a definite x 
fermion occupation number. The Hamiltonian is then 
block diagonal, each block corresponding to a distinct set 
of X fermion occupation numbers. Thus all eigenstates 
in the extended Hilbert space take the factorized form, 

|*> = |Xe;g> ^ \Mg)\g) (9) 
and x\ij)^X(i3)AQ) = '^(-y)J^) (10) 

where rii^ij-^^ = "^'^2°^^ l-^s) is a many body eigen- 
state in the matter sector, corresponding to a given Z2 
field of \Q). In terms of bond fermions, spin operators 
become, 

< = ic^(x{^J)^+X%J)^) (11) 

= (12) 

Three components of a spin operator at a site, gets 
connected to three different Majorana fermions defined 
on the three different bonds ! Written in the above 
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form, the effect of cr° on any eigen-state, which we re- 
fer to as a "spin flip", becomes clear. In addition to 
adding a Majorana fermion at site i, it changes the bond 
fermion number from to 1 and vice versa (equivalently, 
'^(ij")a ~^ ~ '^{ij)a)j cit the bond {ij)a. The end result 
is that one tt flux each is added to two plaquettes that 
are shared by the bond (ij)a (figure 2). We denote this 
symbolically as 



ICi TTi 

{ij)a '^2{ij)a 

(13) 



with TTi(^ij)a and TT2{ij)a defined as operators that add 
TT fluxes to plaquettes 1 and 2 shared by a bond {ij)a 
(figure 2). Further 'I'l^y^a = 1, since adding two tt fluxes 
is equivalent to adding (modulo 2tt) zero flux. 

Now we wish to calculate spin-spin correlation func- 
tions in physical subspace. Since the spin operators are 
gauge invariant, we can compute the correlation in any 
gauge fixed sector and the answer will be the same as in 
the physical gauge invariant subspace. (We have con- 
firmed this by a calculation in the projected physical 
subspace.) So we consider the 2-spin dynamical corre- 
lation functions, in an arbitrary eigen-state of the Kitaev 
Hamiltonian in some fixed gauge field configuration Q, 



sfHt) ^ {Mg\{gwnt)a'mg)\Mc 



(14) 



Here A(t) = e^^^Ae is the Heisenberg representation 
of an operator A, keeping %= 1. As discussed above, 



<^]mQ)\Mg) 

<ytit)\Q)\Mg) 



(15) 
(16) 



where, \Q^°-'^i'^^) denote the states with extra tt fluxes 
added to Q on the two plaquettes adjoining the bond 
{ik)a){{lj)b)) and E is the energy eigenvalue of the eigen- 
state |5|A^g). Since the Z2 fluxes on each plaquette 
is a conserved quantity, it is clear that the correlation 
function in equation ([14]) which is the overlap of the two 
states in equations US! [16]) is zero unless the spins are 
on neighbouring sites. Namely, we have proved that the 
dynamical spin-spin correlation has the form, 

^iji^) — g(ij)A^)^a,b , ij nearest neighbors (17) 



= 



otherwise. 



Computation of 5ij(0) is straight forward in any eigen- 
state \Mg). For the ground state where conserved Z2 
charges are unity at all plaquettes, the equal time 2-spin 
correlation function for the bond (ij)a is given by the 
analytic expression: 



Where cos9{ki, 
louin zone, 



V3 
167r2 



cos^(fci, k2)dkidk2 



BZ 



Ek = yi^^TAf) , in the Bril- 

2(J:rCOsfci -1- Jy C0Sk2 + Jz), Afc = 



2(Ja;Sinfci -I- Jj,sinfc2), ki = k.ni, ^2 = k.n2 and ni.2 = 
^Gj. ± are unit vectors along x and y type bonds. 

At the point, = Jy ^ J^, we get = -0.52. 

To compute g{ij)^{t) we substitute for the cr's from 
equation ([7]) and ([8]). We choose a gauge where U(^ij)^ = 
-1 implying x\ij)^ \Q) = x|ifc), \Q) =0. We note that 
the above conditions imposed at t = will continue to 
be true at all times since the bond fermion numbers are 
conserved. We then have. 



(18) 

The time dependence evolution can be expressed in terms 
of the hamiltonian and noting it is diagonal in the number 
operators x^X: we get, 

g(,,>(0 = (A^e|e^^[e"l*zc.(0)e-^^[e"°l*(~l)c,(0)|Me) 

(19) 

where H[Q^°-] is the tight binding hamiltonian in the 
background of the static gauge field configuration CJ'". 
The (—1) factor is U(ij)a- This expression can be writ- 
ten in terms of the time evolution under H[Q] as follows, 

9{r3)St) = (Me|*C»(t)r(e"2J„/>(,,,„c,(r)c,(r)dr^ 

W(,,)„C,(0)|Me) (20) 

The above equation is written in an arbitrary gauge. 

We have thus derived a simple but exact expression for 
the spatial dependence of the two spin dynamical correla- 
tion function. We have also obtained an exact expression 
for the time dependence in terms of the correlation func- 
tions of non-interacting Majorana fermions in the back- 
ground of static Z2 gauge fields. Equation (j20|l represents 
the propagation of a Majorana fermion in the presence 
of two injected fluxes. It can be treated as an X-ray edge 
problem and computed in terms of the Toeplitz determi- 
nant. We will not do this now but proceed to discuss 
some general features of our results. 




-iHt 



FIG. 2: Time evolution and fractionization of a spin flip at t 
= at site i, into a 7r-flux pair and a propagating Majorana 
fermion. 'Shakeup' of the Majorana fermion vacuum to an 
instantaneous addition at t = 0, of a yr-flux pair is not shown. 



The notion of fractionalization of 
natural interpretation of our results 



spm- 



flip quanta is the 
i, Q . Consider time 
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evolution of a single 'spin-flip' at site i given in equation 
(fTH)) . Using the notation introduced in equation ([13]) we 
have, 

(21 

A spin-Hip at site i at time t = is a sudden perturba- 
tion to the matter (Majorana fermion) sector, as it adds 
two static TT-fiuxes to adjoining plaquettes. The time 
ordered expression represents how a bond perturbation 
term, i2u(^ii.)aJaCiCf. evolves the Majorana fermion state, 
in 'interaction representation'. At long time scale the 
resulting 'shakeup' is simple and represents a rearrange- 
ment (power law type for gapless case) of the Majorana 
fermion vacuum to added static 7r-flux pairs. The Ma- 
jorana fermion, produced by a spin-flip, c,;(t) propagates 
freely, as a function of time. 

As spin-flip at site i is a composite of a Majorana 
fermion and 7r-flux pair (equation 13), two spin corre- 
lation function defines the probability that we will de- 
tect the added composite at site j after a time t. As the 
added tt-Sux pair do not move, the above probability is 
identically zero, unless sites i and j are nearest neighbors 
and spin components are a = b. This is why the spatial 
dependence of two spin correlation functions are sharply 
cut off at nearest neighbor separation. The asymptotic 
response to an added tt-Aux pair and free dynamics of 
the added Majorana fermion control the long time power 
law behaviour of our only non vanishing nearest neighbor 
two spin correlation function. 

Further, for a given pair of nearest neighbor sites, only 
one Ising spin pair of a corresponding component is non- 
zero. Other pairs and cross correlation functions vanish. 
More specifically, for a given bond the only non zero two 
spin correlation function is the bond energy. 

What is unusual is that the above result is true in all 
eigen-states of the Kitaev Model, irrespective of ener- 
gies. It follows that it is valid for thermal averages too. 
This is an unusual result, indicating exact fractionaliza- 
tion occurring at all energy scales. In known models such 
as ID repulsive Hubbard model or spin half Heisenberg 
chain, fractionization is only a low energy asymptotic 
phenomenon. Our results show the all energy scale exact 
confinement of the spin-flip quanta, and exact deconfine- 
ment of the Majorana fermions in the Kitaev model. 

It is interesting to see that the above is a special prop- 
erty of the Kitaev Model. When we perturb it by adding, 
for example, a magnetic field term or make bond terms 
non-Ising, 7r-fluxes acquire dynamics. This means that 
the probability amplitude of finding the composite parti- 
cle intact at a farther site is finite (though exponentially 
small as a function of separation) and not strictly zero. 

Multi spin correlation functions can be calculated in 



our formalism. Further, quantum entanglement, a key 
notion in quantum computation and quantum informa- 
tion, is ultimately connected with some complicated 
multi-spin correlation function. We have computed some 
entanglement measures, but do not discuss them in the 
present paper. 

To summarise, this paper presents certain exact ana- 
lytical results for the spin dynamics and a spin-flip frac- 
tionization scheme for the Kitaev Model. As this non- 
trivial spin model is also a model for topological quan- 
tum computation, our exact results should provide in- 
sights into qubit dynamics and possible ways of generat- 
ing emergent topological qubits. Our formalism, which 
uses the factorized character of the eigen-functions in the 
extended Hilbert space, is easily adapted to the calcula- 
tion of multi-spin correlation functions, which is a key 
step in the calculation and understanding of quantum 
entanglement properties. 
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